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ABSTRACT 



We show how it is possible to give a precise and unambiguous imple- 
mentation of the Witten- Veneziano formula for the rj mass on the 
lattice, which looks like the formal continuum one, if the expression of 
the topological charge density operator, suggested by fermions obeying 
the Ginsparg- Wilson relation, is employed. By using recent numerical 
results from simulations with overlap fermions in 2 (abelian Schwinger 
model) and 4 (QCD) dimensions, one obtains values for the mass of 
the lightest pseudo-scalar flavour-singlet state that agree within errors 
with theoretical expectations and experimental data, respectively. 



1 On leave of absence from Dipartimento di Fisica, Universita di Roma 11 Tor Vergata" , 
Via della Ricerca Scientifica, 1-00133 Roma, Italy. 
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1 Introduction 



Lacking any reliable analytical way to compute bound-states masses in a 
strongly interacting theory, a crucial confirmation of QCD as the theory of 
hadron dynamics has come in recent years from our ability to compute from 
first principles the masses of the lightest hadrons in Lattice QCD (LQCD) |]]. 
This has been done with remarkable success for a number of mesonic and 
baryonic states, both in the quenched approximation (no fermion determi- 
nant) and, more recently, also in the full (unquenched) theory ||. 

A somewhat difficult problem is, computationally, the inclusion of dia- 
grams violating the so-called OZI-rule ||, i.e. of diagrams involving purely 
gluonic intermediate states: these require the calculation of the trace of the 
inverse Dirac operator, and are affected by large statistical fluctuations. The 
very existence of such diagrams is a key property of QCD, making it differ in 
an essential way from any naive quark model. Since OZI-violating diagrams 
are the ones that distinguish flavour-singlet from flavour-non-singlet states, 
their inclusion in LQCD is crucial and should allow for the theoretical com- 
putation of subtle phenomena, such as p — uo — <f) and tt — r\ — rf splitting 
and mixing. In fact, one of the most striking experimental facts in hadronic 
physics is the contrast between the vector (as well as the tensor) mesons, 
that appear to satisfy the OZI-rule to high accuracy, and the pseudo-scalar 
mesons, whose (mass) 2 eigenstates are very far from being ideally mixed 
(the OZI-prediction), and whose eigenvalues exhibit a peculiar hierarchical 
structure, m 2 <C m 2 «C m 2 ;. 

The technical problems mentioned above have so far prevented in LQCD 
a fully-fledged evaluation of the r/ mass at the same level of accuracy as that 
reached for the other light hadrons, although noticeable progress has been 
recently made in this direction |4], |5j . 

In quenched LQCD the first evaluation of the r/ mass Q dates back to the 
work of ref. 0, where the OZI-violating (ZV) and the OZ I- conserving (ZC) 

2 Throughout this paper we will indicate by rj the flavour-singlet pseudo-scalar meson, 
ignoring vector flavour symmetry-breaking effects. We will work in Euclidean metric with 
hermitian 7 matrices satisfying {7^,7,,} = 26^ and tr (75 71 72 73 74) = 4ei234 = 4. Left- 
handed and right-handed fermions are defined through the equations 75 ipL — and 
Ibi'R = — ^Ri respectively. Colour generators, T a (a = 1, ...,N% — 1), in the fundamental 
representation are normalized according to tv(T a T b ) = S ab /2. 
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contributions to the 2-point meson-meson Green function 

r 5 M = /df(p°(f, t) p»(o)>| q _ M (i) 



were computed. In eq. 



N f 

P° = Y, ^5</V (2) 



r=l 



is the pseudo-scalar singlet quark density, with Nf the number of light quarks. 
The (quenched value of the) rj' mass is identified as the amount by which the 
singlet pseudo-scalar gg-boson pole gets shifted owing to the iteration of 
Nambu-Goldstone (NG) particle exchanges, the sum of which builds up the 
full (unquenched) expression of F. This leads to the formula 

r zv (t) 

<' = 2m *}™ jtjrf(*) • ( 3 ) 

Numerical data were encouraging, but far from convincing. Although con- 
siderable progress has been achieved along these lines in more recent simu- 
lations 0, ||, in our opinion one cannot regard the present numerical results 
as really conclusive. 

From a more theoretical point of view, leaving aside the lattice for a 
moment, the striking difference between vector and pseudo-scalar mesons 
appeared as a serious problem within the generally accepted framework ac- 
cording to which the light pseudo-scalar mesons are the (quasi) NG bosons of 
a spontaneously (and also explicitly) broken chiral symmetry. As first noted 
by Glashow ||, a naive use of spontaneously-broken chiral symmetry calls 
for nine - rather than eight - light NG bosons (Nf = 3). Thus one would 
naively expect that chiral symmetry should "protect" the OZI-rule in the 
pseudo-scalar channel, while the opposite is experimentally true. 

Awareness of the problem was reinforced by a very influential paper by 



Weinberg |L0j in which, under what he calls "plausible assumptions", the 
bound 

^octet,/=0 < 

V3m w (4) 

was derived in the absence of {7a (1) anomaly contributions. 

Since the very beginning, the famous ABJ anomaly (Tl| was suspected 



to have much to do with the resolution of the {7a(1) puzzle, as the strong 
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(QCD) axial anomaly appears only in the flavour-singlet channel. But the 
fact that the anomaly itself is a total divergence (albeit of a non-gauge- 
invariant current) made it look irrelevant for the solution of the problem. 



Kogut and Susskind |12| were the first to notice that the ABJ anomaly could 
be made effective again, if one assumed the existence of a "ghost" particle, 
a massless state coupled to the non-gauge-invariant current, but decoupled 
from any gauge-invariant operator. 

Things took a decisive turn when, in 1976, 't Hooft |L3j pointed out that 
the resolution of the Ua(1) problem had to be related to the existence of 
topologically non-trivial gauge field configurations (instantons) in Euclidean 
QCD. Soon after, Crewther emphasized that 't Hooft's crucial observa- 
tion could only be considered a solution of the problem, if it could pass a 
number of theoretical tests, in particular that of satisfying both the con- 
ventional (i.e. non-anomalous) and anomalous Ward-Takahashi Identities 
(WTI's). Several puzzles, related to the dependence on the number and 
masses of quark flavours appeared, in fact, to be still present fl4j . 

Fortunately, these questions can be addressed and solved by exploiting 
the simplifications that occur in certain limits of QCD. One such limit, em- 



ployed in ref. [15 



is the 't Hooft limit (N c — > oo, with g 2 N c and Nf held 
fixed [16|]). The other possibility, proposed in ref. |L7||, proceeds by assuming 
that anomalous flavour-singlet axial WTIs retain their validity order by order 
in an expansion in u = Nf / N c around u = 0. In either cases one can derive 
the leading-order Witten-Veneziano (WV) relation 



2N 



A. 



(5) 



where F n ~ 94MeV is the pion decay constant and A is the "topological 
susceptibility" . A is formally defined by the equation 



A = j d 4 x (Q(z)Q(O)) 
where Q{x) is the topological charge density 



YM 



(6) 



Q(x) 



9' 



32tt 



2 ^fivpa 



tr[F^F pa (x)} 



(7) 



3 This expansion is not to be confused with the topological expansion (N c — > oo at fixed 
g 2 N c and u), introduced in ref. Jig], 
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and the trace is over colour indices. The notation (...)|ym in eq. (|6|) means 
that the QQ-correlation function is to be computed in the pure Yang-Mills 
(YM) theory, i.e. in the absence of quarks. 

We recall that, within the approach of ref. |15| , eq. (Q) is expected to 



be valid for large N c , while, relying on the derivation given in ref. [fLTA , one 
would conclude that the rj' mass formula is valid for any value of N c at leading 
order in u. We also notice that the limit u — > at fixed g 2 N c encompasses 
the 't Hooft limit (N c — > oo at fixed Nf) and that the equivalence of eqs. (0) 
and (|^) essentially relies on the validity of anomalous flavour-singlet WTIs 
order by order in « (17]. 



Since the pioneering works of refs. |L3|, and [17]], a large number of 



papers (see, for instance, |r9|-[]25| and references therein) have appeared, 
aimed at checking the validity of eq. (]5|) by comparing the number obtained 
by inserting in it the value of A extracted from pure gauge lattice simulations 
with the experimental magnitude of the rj mass Q. 

The formal eqs. (J5|) and (D), when translated in any regularized version 
of QCD, such as LQCD to which we will refer in the following, become more 
complicated and quite a number of subtleties have to be dealt with in order 
to be able to determine the correct field theoretical expression of A, which 
should be inserted in the rf mass formula. 

Two main problems need to be solved to make eqs. (|5|) and (|6|) rigorous 
and of practical use. One has to i) find a properly normalized lattice defini- 
tion of the topological charge density Ql] ii) subtract from the Ql(x)Ql{0) 
operator product appropriate contact terms, as required to make it an inte- 
grable (operator-valued) distribution. 

Three types of prescriptions have been proposed upon addressing the first 
problem. They can be succinctly described as follows. 

• Take any naive lattice discretization of F^ v |T{| [|27| |^8[ and construct 
Ql through (the lattice version of) eq. (0) by properly normalizing 
it& ' 

• Define Ql by making reference to the anomalous flavour-singlet axial 
WTIs p9| M M- 



4 The question of how eq. (||) should be modified beyond the chiral limit was addressed 
in refs. @, 
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Exploit the topological nature of Q to give a definition of it appropriate 



to a discrete manifold 31 



The second problem is actually more subtle than stated above. In fact, it 
is not sufficient to subtract (possible) divergent contact terms from the r.h.s. 
of eq. (|3|), by taking properly into account the mixing of Q(x)Q(0) with lower 
dimensional operators. It is also necessary to fix finite contact terms of the 
form, say, c8(x), since, as also stressed in ref. [p4| , they would contribute a 
constant to the if mass when inserted in eq. (|6|). 

The question of the subtraction of lower dimensional operators was ad- 
dressed in an incomplete way and in perturbation theory in ref. ||19|| , and in 
a more refined, non-perturbative, way in ref. |[22|| . The issue was taken up 



again more recently within the overlap formulation of LQCD |32|, |30[ (see 
ref. |33[ for a recent review). 



The purpose of this paper is to give a rigorous derivation of a formula for 
the rj mass, which could be unambiguously used in numerical simulations. 
The validity of such a formula is a crucial test of QCD, as it directly relates 
the non-vanishing of the rj mass in the chiral limit to the explicit breaking 
of the Ua(X) symmetry induced at the quantum level by the gluon anomaly. 

Our strategy is to use (lattice) regularized anomalous flavour-singlet axial 
WTIs of full QCD in order to properly construct the (QQ) correlation func- 
tion, as a well defined (integrable) distribution. The remarkable result of our 
analysis is that no subtraction is needed to this end in the chiral limit, if the 
definition of lattice topological charge density suggested by fermions obeying 
the Ginsparg- Wilson (GW) relation [[34j] (e.g. overlap fermions [3(J) is 



employed. In this formula for the rj' mass, which looks like the naive 

continuum one (i.e. which closely parallels the structure of eqs. (|) and (^|)), 
can be obtained. 

The case of Wilson fermions is somewhat more involved and will be dis- 
cussed in a forthcoming publication f35 |. 



The paper is organized as follows. In sect. |2| we recall the derivation 
of the rj mass formula in the continuum. In sect. ^ we show how one can 
give a precise meaning to this relation in the QCD regularization offered by 
GW fermions. In sect. |] we discuss the numerical evidence in favour of the 
WV formula in 2 and 4 dimensions, which emerges from the existing overlap 
fermion simulation data. Conclusions can be found in sect. [5[ In an Ap- 
pendix we derive several relations between unquenched and quenched QCD 
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quantities, involving moments of the topological susceptibility, including the 
WV formula. 

2 The rf mass formula in the continuum 

In this section we review the formal derivation of the WV formula in the 
continuum. We start from the anomalous flavour-singlet WTI 

d,(Al(x)Q(0)) = 2m q (P°(x)Q{0)} + 2N f (Q(x)Q(0)) , (8) 

where 

N f 

K = E Vi^r (9) 

r=l 

is the axial current, P° was defined in eq. (§), Q is the topological charge 
density and m q is the quark mass. It is understood that all operators in 
eq. (§D are finite (i.e. have finite insertions with any string of renormalized 
fundamental fields), so that every term in eq. (|9]) is finite as soon as x ^ 0. 
The Fourier transform of eq. (§) in the chiral limit reads 

ip^J d 4 xe-^(A°(x)Q(0)) + CT(p) = (10) 
= 2N f J d 4 xe- ipx (Q(x)Q(0)) + CT(p) = 2N fX t(p) ■ 

A contact term, CT(p), has been added to both sides of eq. ( |10|) to make 
them separately finite. CT(p) is a polynomial of degree 4 in p, which vanishes 
at p = 0, because the Fourier transform of the l.h.s. of eq. (H) is certainly 
finite (actually zero) at p = in the full theory. As will become clear in 
the following, CT(p) plays no role in the main line of arguments we present 
in this paper, since finally we will be only interested in the value of Xt at 
p = O0. 

The absence of zero-mass particles in the pseudo-scalar singlet channel 
implies the "sum rule" 

Xt(0) = 0. (11) 

5 The presence of CT(p) will, instead, affect the determination of quantities like Xt(0) = 
dxt(p 2 )/dp 2 \ p =o, unless CT(ff) oc p 4 Q. Xj(0) is relevant to the discussion of the so-called 
proton spin crisis problem ||37| . We will elaborate further on this issue in the Appendix. 
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This relation is the famous statement that the topological susceptibility in 
full QCD vanishes in the chiral limit. 

One can derive from it a formula for the rj' mass whenever the topological 
susceptibility in the absence of quarks is non-zero. In this situation, in fact, 
the vanishing of Xt{p) at p = in full QCD to any order in u, can possibly 
take place only if there exists a qq- meson contributing to Xt(p), whose mass 
goes to zero (linearly) as u — ► |15| |17|| . Under these circumstances the 
limit p — > and the expansion of Xt{p) around u = do not commute. In 
the Appendix we give a derivation of the WV formula along these lines and 
prove a number of other interesting related equations. The formulae that 
are obtained are derived under the assumption that at p ^ all the relevant 
quantities and correlators are smooth functions of u and, in particular, that 
taking the limit u — > is equivalent to neglecting the fermion determinant. 
In the following, for short, we will refer to this assumption as the "smooth- 
quenching hypothesis" . 

To leading order in u one then obtains for the rf mass either the formula 



^/^a-WW))! (12) 



o AT f 

ml = -jf J d*x (Q(x)Q(0))\ ym . (13) 



or, equivalently 



We stress that, under the above assumptions, the r.h.s. of ([13]) (as well as 
that of eq. (|12|)) is ultraviolet-finite, because CT(p) is a polynomial in p that 
vanishes at p = for any value of u. 

In the next section we will show how it is possible to give a precise meaning 
to the above equations, starting with the chirally invariant regularization of 
QCD offered by fermions obeying the GW relation. In this framework we 
will be able to derive an unambiguous formula for the rf mass, immediately 
usable in actual Monte Carlo simulations. 

We conclude by remarking that eq. fli"2] ) cannot be employed for direct 
numerical evaluations, because its r.h.s. is actually identically zero at finite 
volume, as there are no massless particles in a finite box. However, it can be 
put in a form that is amenable to numerical simulations, if one goes one step 
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further and rewrites it in the form f] 



2N f 



<o\QW) 



quenched 



(14) 



This equation follows from the fact that in the quenched limit the integral in 
the r.h.s. of eq. (|T2|) is exactly given at the chiral point by the contribution 
of the pole of the singlet pseudo-scalar particle, as the mass of the latter 
vanishes in that limit. No explicit reference to the size of the volume in 
which the system is enclosed appears anymore in eq. (|14|) . We will make use 



of these considerations in sect. 3.1. 



3 Ginsparg- Wilson fermions 



Regularizing the fermionic part of the QCD action using GW fermions 34, 32 



3T| offers the great advantage over the standard Wilson discretization ]I[ that 



global chiral transformations can be defined, which are an exact symmetry of 
the massless theory, as in the formal continuum limit. This is a consequence 
of the relation @ [] 

l5 D + D l5 = D l5 D , (15) 

obeyed by the GW Dirac operator, D. In this regularization the Ua(1) 
anomaly is recovered a la Fujikawa []3~8[| , as a consequence of the non-trivial 
Jacobian that accompanies the change of fermionic integration variables in- 



duced by a Ua(X) lattice transformation. In fact, it can be shown |39|] that 
under the infinitesimal global transformations 

1> = ij/ + ery 5 (l - D)if/ 

^ = ^ + e^ T 5 (16) 



the massless action, thanks to (|15T) , remains invariant, while the functional 
integration measure gets multiplied by the factor 

J = exp{e J d A x Ti\^D{x, x)]} , (17) 



6 See the Appendix for the relevant normalizations and notations. 
7 We set the lattice spacing, a, equal to 1. It is, however, understood that in all lattice 
formulae the continuum limit should be finally taken. 
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where Tr is over colour and spin indices. In view of these results it is con- 
cluded that 

• the index theorem for the GW Dirac operator reads (for each flavour) 
n R -n L = index(D) = - J d*x Ti[y 5 D(x, x)} ; (18) 



for sufficiently smooth gauge configurations |42] one has 



^Tr[ l5 D(x,x)] = ^ w tr[F^F pCT (x)] ; (19) 
Z 6271 



• the anomalous flavour-singlet WTIs in the presence of Nf massless 
fermions take the form 

O AT r- 

= ^Jd i x (Tr[ l5 D(x, x)}6) + (6 A 6) , (20) 

where O is any finite (multi) local operator. In eq. (p0| ) we have sym- 
bolically indicated by e8 A the infinitesimal variation of O under the 
transformations (|T6|). Equation (^) assumes the absence of a £7a(1) 
NG boson. 

Starting with the local version of the transformations (|T(|), one gets in 
the chiral limit the local WTIs 



O AT 

VM»6) = ^(Tr[ 75J D(x, x))6) + (5 A 6) , (21) 



where 5 A is the local variation of 0\ V M v4°(x) is the divergence of the 
singlet axial current, a quantity that formally vanishes when integrated over 
all space-time. 

Neither Tr [y 5 D(x, x)] nor A ^ (x) are finite operators, but finite linear com- 



binations are easily constructed [H3|, PHI]. In fact, since the second term in 



the r.h.s. of eq. (p20|) is obviously finite, if O — > A n = vp R (xi)...?pji(x n ), it 
follows that the integrated operator / d A x Tr [75.0(0;, x)] is also finite, as it 



8 Actually the GW relation (15) alone may not be enough to ensure the validity of 
eqs. ( |l8|) and ( [l9| ) below | f40| . These equations are, however, valid if overlap fermions are 
used pl[ and, more generally, if the admissibility condition of ref. G2] is fulfilled. 
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has finite insertions with any string of renormalized fundamental fields. This 
statement in turn implies that Tr[75.D(x, x)] can only mix with operators of 
dimension < 4 and vanishing integral, hence only with V^A^x). Calling Z 

the relevant mixing coefficient, one can define the finite operators Q and A ^ 
by writing 

Q(x) = \tt[ 15 D(x, x)} - ^»A%x) (22) 

A^x) = (1 - Z)Al(x) . (23) 

Z is logarithmically divergent to lowest order in perturbation theory and 
vanishes as u — > 0. Its perturbative expansion starts at order g 4 . 

With the above definitions the renormalized singlet axial WTI becomes 

V,(A°Jx)6) = 2N f (Q(x)6) + {5\6) . (24) 



We expressly note that 1) as the comparison of the WTIs (£l]) and ( p24|) shows, 
the operators Q and A° are already correctly normalized, so no further (finite) 
multiplicative renormalization is required; 2) the line of arguments given 
above also proves that there is no mixing of Tr[7 5 D(a;, x)] with operators of 
dimension smaller than 4 (like the pseudo-scalar quark density), which would 
bring in dangerous power divergent mixing coefficients. This fact is related 
to the absence of power-like divergent, 1/a, subtraction in the definition 
of renormalized quark mass. All these are very distinctive features of GW 
fermions which are at variance with what happens in the case of Wilson 
fermions |. 



3.1 The rj mass formula on the lattice 

A formula for the rj mass can be easily obtained from the WTI (|21|), assuming 
that in the limit u — > the rj' mass is 0(u). A more sophisticated derivation 
can be found in the Appendix. 

One starts by defining the lattice Green function 

XtM = 2Ar/ d 4 *e-^V M (iJ(*)g(0)> , (25) 

9 We are indebted to M. Liischer for pointing out to us the possible existence of mixing 
between Q and P° in the case of lattice Wilson fermions. 
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with Q given by eq. (|22|). The same observation, that was made in sect. |2] 
after eq. (|10D about the presence of counter-terms necessary to make the 
correlator (^) finite at p 7^ 0, applies also here. However, as this counter- 
term does not play any role in the determination of the 7/ mass, for brevity 
we will not indicate it explicitly in the formulae of this section. 
In the full theory, where the 77' is massive, one gets 

Xti(0) = 0, (26) 

which is the (lattice) regularized version of eq. (|Tl|). 

The key observation at this point is that, at vanishing quark mass, in 
the limit u — > 0, only the 77' pole contributes to Xth(p), as explained below 
eq. (JI3]), giving (see eqs. ( pE9"| ) and ( J5Bp of the Appendix) 

Fl 



lim lim XtM 

p— >0 u— +0 



u=0 



(27) 



Recalling that F n = 0(vAQ [16], to leading order in u, we see that 
consistently the combination in the r.h.s. of this equation has a finite limit 
as u — > 0. 

On the other hand, using the WTI (|2ll) with O = Q, we can put eq. (g^) 
in the form 



X*l(p) 



(28) 



Z) J rf 4 xe-^(^Tr[ 75J D(x,x)](iTr[ 75J D(0,0)]-^V,4(0))) 



since = 0. Taking, as before, the limit p 
we get 

F 2 

T 

— m„/ 

M=0 



iriZ 

2N- n 



lim lim 

p-»0 



after u -»• in eq. (EH), 



Z) / rf 4 xe-^(-Tr[ 75 .D(a;,a;)]-Tr[7 5 D(0,0)]) + 



1 



m„, 

Z2N t " 



(29) 



where eq. (^7|) has been used in the l.h.s. and the last term comes from 
the contribution of the rj' pole to the second term in the r.h.s. of eq. 
Recalling that Z vanishes when u — > 0, we obtain 



u=0 p^O u— >0 



1, 



lim lim / ci 4 a;e- 4pz (-Tr[75 J D(a;,a;)]-Tr[75 J D(0,0)]) . (30) 
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Under the smooth-quenching hypothesis the limits u — > and p — > in this 
expression can be readily performed in the order indicated, as the first one 
simply amounts to setting the fermion determinant equal to unity. One gets 
in this way 



Fl 2 -/d 4 x(^IV[ 75J D(x,x)]ilV[ 75J D(0,0)]>| YM . (31) 



2N f -* 



u=0 



The restriction to the pure YM theory, indicated in eq. fl3~T|), is an obvi- 
ous consequence of the fact that, for a Green function with the insertion of 
only gluonic operators, neglecting the fermion determinant is tantamount to 
limiting the functional integral to the pure gauge sector of the theory. 

A relevant question at this point is to ask for which value of N c eq. fl3T|) 
is supposed to be valid. The answer depends on the detailed behaviour of 
QCD with Nf. Various possibilities are logically envisageable. 

1) In the most favorable situation, in which the limit u — > of Xti{p) 
exists and is equal to the value it takes at Nf = (fermion determinant 
equal to 1), the formula fl3~T|) is valid for any value of N c and for each N c it 
yields the mass of the rj' meson to leading order in u in the world with the 
corresponding number of colours, without any 0(l/iV c ) corrections [17 . 



2) If quenching can be attained only in the limit in which the number of 
colours goes to infinity, then eq. (|3~I| ) will yield a formula for the rriii valid 
up to 0(l/iV c ) corrections [|15| . 



3) Finally it may happen that taking u — > does not correspond to 
quenching. In this case one cannot pass from the fairly complicated eq. ( p0|) 
to the more useful formula (j3T|) . 

Equation fl3T|) has precisely the same structure as the naive formula (|5|), 
when the expression (|j) of the YM topological susceptibility is inserted in it. 
Remembering the meaning of | / d 4 xTi[y 5 D(x, x)] (eq. (|I8|)), eq. (|5lD can be 
rewritten in the suggestive form 



2 



VTil 



2N f " 



= lim ^-"^)> , (32) 



where ((ur — til) 2 ) is the expectation value of the square of the index of the 
GW fermion operator, D, and V is the physical volume of the lattice. 

It may be interesting to see what are the implications of the hypothetical 
situation in which the rj' mass would not vanish (linearly) with u. In this 
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case the two limits u — > and p — > of xtr, would commute, leading to the 
chain of relations 

= lim lim Xth{p) = hm lim Xth{p) = lim ^ Ur l ^ - > , (33) 

«->0p^0 p^O V— *oo v 

which imply either that the non-negative quantity (tir — ul) 2 vanishes iden- 
tically for any gauge configuration or, at least, that the quantity {{n R — n L ) 2 ) 
grows less than linearly with the number of lattice points. 



4 Numerical results 

In the form (|32|) the r( mass formula can be directly compared with overlap 
fermion simulation data, because the problem is reduced to counting the 
number of zero modes of the overlap-Dirac operator and measuring their 
chirality. There exist data both in 2 (abelian Schwinger model |]45|| ) and in 
4 (SU(3) gauge theory |I(J) dimensions. 



4.1 2 dimensions: abelian Schwinger model 



The validity of the WV formula can be explicitly checked in the case of 
the two-dimensional abelian Schwinger model |47| , since the model is exactly 
soluble |i7, i5fl. Numerical checks are also possible, as there exist very recent 
simulations that make use of the overlap fermion formalism 4^ . 



For the quantities of interest, one finds the formulae ES 



4:71 



Xt(p) 



47T 



e 2 /7r 



p- 



! /tt 



(34) 



(35) 



These formulae are for Nf = 1. The generalization to the case with Nf 
species of fermions can be found at the end of the Appendix. 

By inspection one immediately recognizes that Xt(p — 0) = and that 
there exists one single particle in the spectrum, a boson (the analogue of the 
rj in QCD) with mass 



71 



(36) 
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The residue of the rj' pole is given by |49 



P 2 



mw)\ 2 = (-^=f . 

v 7TV47T / 



(37) 



As the negative term in eq. (|35|) is fully contributed by the fermion determi- 
nant, the rj' mass formula in 2 dimensions reads 



2 

mi, 



Xt(0) 



quenched 47T 2 



(38) 



Comparing eq. fl38|) with (^) , one immediately checks the consistency of (|38|) 
with (|36"l). Incidentally, looking back at eq. (j31~l), we conclude that = 

1/V27T- 

Numerically one finds ((tir — Ul) 2 ) = 2.7 ± 0.4, at 1/e 2 = 6 in a lattice 
with iV = 24 2 points [ J5| , implying 



(K-n L ) 2 ) 2.7 ±0.4 



(4.7 ±0.7) x 10~ 3 . 



(39) 



N 24 2 

Although this number compares fairly well with the expected result (eq. (|38|) ) 



1/6 



471-2 4^2 



4.2 x 10 



-3 



(40) 



further simulations at different /3's and larger volumes are needed to have 
properly under control all sources of systematic errors. 



4.2 4 dimensions: SU(3) gauge theory 

Using the spectral flow method fl50|, the expectation value of {n R — n L ) 2 in 
an SU(3) gauge theory has been computed for three lattice volumes (N = 
8 3 x 12, 8 3 x 16, 16 3 x 8) at /3 = 5.85 in the first paper of ref. gf|. The 
data show a fairly good scaling with N. If the value of the string tension 
a = (440 ± 38MeV) 2 |2| is used together with = 0.2874(7) || to fix 
the lattice spacing, one finds (V = Na 4 ") 

A = ^ R ~ Ul)2) = (198 ± 20 MeV) 4 , (3 = 5.85 . (41) 
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The same authors 
volumes. At (3 
get 

((n R -n L f) 



have also produced data at larger values of f3 and 
6.0 and iV = 16 3 x 32, using a^/a = 0.2189(9) [0, they 



.4 



V 



(l94±20MeV) 4 , (3 = 6.0 



(42) 



Errors in eqs. ( fil]) and fl42l) have been computed by combining in quadrature 
the uncertainties on ((tor — tol) 2 ), <t and ax/a. It is remarkable that the 
central value of A is quite stable with f3 and very near to the number required 
to match the actual value of the rf meson mass [pjj . Clearly more work is 



needed (larger (3 values and volumes) before these encouraging indications 
are fully confirmed and the magnitude of the error is properly assessed. 



5 Conclusions 

In this paper we have given a derivation of the WV formula for the rf mass 
in lattice QCD, by exploiting the flavour-singlet axial WTIs of the theory. 
If fermions obeying the GW relation ( |T5D are used, there exists a natural 
definition of topological density (eq. (|1^)), for which the naive form of the 
WV formula holds without the need of introducing any subtraction. 

The validity of a formula of the WV type for the rf mass, i.e. of a formula 
that is telling us that the 0(A QCD ) contribution to originates from the 
breaking of the Ua(1) symmetry due to the gluon anomaly, is a key test of 
our understanding of strong interaction dynamics, as described by QCD. It is 
reassuring that such a formula holds even beyond the formal continuum- like 
derivations given in refs. [1^, [I7| . 

In the literature essentially two kinds of approaches have been proposed 
to deal with the problem of computing of the topological susceptibility on the 
lattice, which have led to numerical values for A as good as those obtained 



here. The first one is based on a direct field theoretical definition of A p2 
that takes into account the need for the renormalization of Ql and the sub- 
traction of the operators F 2 and 1 in the short distance expansion of QlQl- 



The second one makes use of the notion of "cooling" [52] to carry out the 
necessary operations of renormalization and subtraction. Both methods are 
reliable to the extent that they are able to capture the topological properties 
of the gauge field configurations that determine the number of zero modes 
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of the Dirac operator. Simulations based on the geometrical definition of 
ref. |H| have not yet led to comparably good results |21 . 



We conclude by noting that a completely analogous set of arguments, 
as those developed in this paper, could be carried out for the pseudo-scalar 
meson (sort of rj/n particle) belonging to the lightest supermultiplet of the 
M = 1 Super- YM gauge theory, the only difference being that in this case 
one would be dealing with Weyl fermions in the adjoint representation of the 
gauge group. In order to assess the numerical importance of gluino loops 
for the restoration of supersymmetry, it would be interesting to compare the 
values of the masses of the other two partners of the lowest lying supermulti- 
plet, as they are obtained in quenched simulations |53], with the number one 



would derive for the mass of the ri/ir particle following the present analysis. 
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Appendix 

We present here a simple (and to our knowledge novel) derivation of a num- 
ber of relations between certain hadronic observables and various quantities 
related to (the continuum limit of) the "topological susceptibility correlator" 
Xt(p) (see eq. (g|)). 

The dispersion relation for Xt(p) 

As Xt(p) is a quantity of mass dimension 4, it will satisfy a three-times- 
subtracted dispersion relation Q that we write in the form 

Xt(p) = Xt(0)+ X 't(0)p 2 + lx"mP 2 ) 2 + (p 2 ) 3 I, (43) 



10 The same results as those obtained below follow, also if one starts with a dispersive 
integral having a smaller number of subtractions. 
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where for a properly defined correlator we have (see eq. (|26|)) 



Xt(0) = 0. 



(44) 



In eq. (H3) / represents the dispersive integral involving lmxt(p)- It is useful 
to express / by separating out the rf contribution from the rest. We thus 
put 



1 



(mi 



) 3 p 2 + m 2 , 



+ 1 



(45) 



where — Bl, is the full physical residue of the rj' pole, which, as explicitly 
indicated, is negative-definite in Euclidean metric. This is a consequence of 
the fact that, by analytical continuation from the Minkowski region jl9], the 
whole dispersive integral is negative-definite. This circumstance solves the 
famous sign problem discussed in |T5 . 



The key observation at this point is that for the vj contribution the limits 
p — > and it — > do not commute, since (if) m 2 , = O(w). We then proceed 
by computing the two sides of eq. fl43[ ) by taking u — > before p — > 0. 

In doing so, in the l.h.s. of eq. (]43| ) we simply get the quantity Xtip) m 
the quenched theory, i.e. 



Xt 



quenched 



limvi(p) . 

u— >0 



(46) 



This follows from the assumption that, since p ^ 0, taking the limit u —>■ 
(before p — > 0) is equivalent to neglecting the contribution of the fermion 
determinant (smooth-quenching hypothesis). In the r.h.s. we expand at 
fixed p 2 the rj' contribution in powers of m 2 ,/p 2 = 0(u/p 2 ) and then we match 
the resulting p 2 power behaviour with the Taylor expansion of Xt{p)\ quenched 
around p 2 = 0. We get in this way the relations 



X"(0)| 



quenched 



lim 



x"(0) + 2-^' 



mt 



(47) 



Xt(0) | quenched 

Xt(0)| quenched 



lim 



lim 



Z?2 

xKo) - 



m 



Xt(0) + 



m\, 



(48) 
(49) 
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Since J is a smooth function of u and p, which vanishes as p — > with its first 
two derivatives, it does not give any contribution to the above equations. 

The smooth-quenching hypothesis implies that the combinations appear- 
ing in the r.h.s. of eqs. (f47|) to fl49|) have a finite limit when u — > 0. As a 
consequence, they can be rewritten in the more expressive form 

R 2 

= X*(0)|quenched ~ *t(0) + 0(u) , (50) 

~ = XiWUnchcd - X' t (0) + 0(u) , (51) 

2^ = Xt (0)| q uenched " x'M + O(u) . (52) 

In this way it clearly appears that, to leading order in u, the physical quan- 
tities in the l.h.s. of these equations are not affected by the presence of any 
possible counter-term, introduced to make the topological susceptibility cor- 
relator finite at p ^ (see the discussion after eq. fllDp). This immediately fol- 
lows from the observation that the w-independent part of such counter-term, 
even if non zero, will cancel between quenched and unquenched quantities in 
the r.h.s. of eqs. (|50D, 



51) and 



The T]' mass formula 



Owing to eq. ([BP, eq. (|50[), taken in the limit u — > 0, is nothing but the WV 
relation. To see this we recall that the residue of the rj' pole is given by the 
formula 



(0\QW) . (53) 

To get an explicit expression of R 2 , , it is sufficient to write the matrix element 
of Q appearing in eq. (^) in terms of the corresponding matrix element of the 
divergence of the flavour-singlet axial current, using the anomalous PCAC 
relation. As we are in the chiral limit, we get 

2N f (0\Q\rf) = <0|VM>'>. (54) 



We can express the r.h.s. of eq. ( |54"D in terms of F n , since, for the sake of 
deriving the WV formula, we are actually only interested in computing R 2 , 
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in the limit u — > 0, where the pseudo-scalar meson decay constants, F w and 
Frf, coincide. Combining the equation []] 



which defines F v > } with 



<0|VM>'> = pN^ml, 
we obtain 



rat 



u=0 



2N f 



u=0 



2N f 



u=0 



(55) 



(56) 



Recalling that F n = O(v^Vc) 0, |18| to leading order in u, we see that 
we consistently have R^/rn^, \ U=Q = 0(1). If we now put together eqs. (|50|) 
and (0), we finally get 



2N 

m v' = -^"Xt (0)| quenched + 0(w 2 ) . 



(57) 



It is amusing to notice that eliminating R^, / m^, between eqs. ( pi]) and (|52l) 
provides a new formula for the r( mass, namely 



Xt 



quenched 



quenched 



xm 



(58) 



It 



which, as indicated above, is more accurate than the WV formula ( [3 
should be stressed that from eq. ( |58| ) m 2 , is expressed, up to 0(u 3 ), in terms 
of differences between "quenched" and "full" quantities that are not affected 
by the presence of possible counter-terms in the definition of the topological 
susceptibility correlator at p ^ 0. 

n The numerical factor y/2Nj in equation ( p5| ) needs some explanation. The factor \/2 
comes from the fact that F„ is defined through the equation (0| VM^° \tt°) = i^m^, where 
tt° = 2- 1 / 2 (uu - dd), Af = qj^X^q and V" = A 3 /2, with A 3 the usual Gell-Mann 
matrix. The factor y/Nf is due to the fact that the flavour-singlet axial current, A ^, is 
a sum over all light flavours with weight 1, while the wave function of the flavour-singlet 
meson has a 1/ \/Nf factor when expressed in terms of qq states. With these definitions 
is of 0(1) in Nf to leading order in N c . 
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A formula for x£(0) 



An interesting formula for Xt(®) can b e derived from eq. (pTJ). In fact, recall- 
ing eq. fl52|), the former can be identically rewritten in the more expressive 
fashion 



2N, 



+ 



rati 



2N, 



+ X' t (0)\ 



quenched 



+ 0(u) 



(59) 



where the first term is 0(l/u) and the term in square brackets is of 0(1), 



as follows from eq. (56). An expansion of this kind for Xt(0) can be useful 
in connection with the discussion of the proton spin crisis problem, given 
in ref. [37], where it was noticed that, in the limit in which Xt(0) is iden- 
tified with F%/2Nf, one recovers the Ellis- Jaffe sum rule [54}|. It would be 



interesting to see whether the correction in square brackets goes in the right 
direction for solving the proton spin puzzle. 

For the study of the proton spin problem one is directly interested in x't{®) 
itself and not in "quenched" -"full" differences, like the combinations fl50|), 

for the rf mass. 



fl5"T|) and fl5"2|) appearing in the formulae (5TD and 
Thus the question arises whether the quantity Xt(0) can be defined in a 
way amenable to unambiguous numerical simulations. This is the case only 
if terms proportional to p 2 are absent in the polynomial expansion of the 
counter-term, CT(p). However, in a (lattice) regularized theory the presence 
of a quadratically divergent term of the type p 2 l/a 2 cannot be excluded and 
will thus have to be compensated for. This means that there is no obvious 
way to define Xt(0) m a wa Y suitable for numerical simulations, at least from 
these simple considerations only. 



The two-dimensional case 

We conclude by noticing that in the two-dimensional abelian (iV c = 1) 
Schwinger model with Nf species of fermions, one can directly verify the 
validity of eqs. (^) to (0), which, in fact, hold exactly without 0(u = Nf) 
corrections. It is sufficient to observe that the generalization of eq. (|3f|) to 
arbitrary Nf reads 



while the ratio (B3) is independent of N 



/■ 
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